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Qh! The one-loop effective potential for gauge models in static de Sitter space at 

finite temperatures is computed by means of the (^-function method. We found a 
simple relation which links the effective potentials of gauge and scalar fields at all 
temperatures. 

In the de Sitter invariant and zero-temperature states the potential for the scalar 
electrodynamics is explicitly obtained, and its properties in these two vacua are 
compared. In this theory the two states are shown to behave similarly in the regimes 
of very large and very small radii a of the background space. For the gauge symmetry 
broken in the flat limit (a — > oo) there is a critical value of a for which the symmetry 
is restored in both quantum states. 

Moreover, the phase transitions which occur at large or at small a are of the first or 
of the second order, respectively, regardless the vacuum considered. The analytical 
and numerical analysis of the critical parameters of the above theory is performed. 
We also established a class of models for which the kind of phase transition occurring 
depends on the choice of the vacuum. 
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1 Introduction 



The aim of this paper is to investigate the properties of the finite-temperature (FT) 
effective potential of the gauge models in the de Sitter space. The temperature [3~ l is 
introduced as the temperature of quantum fields which are in thermal equilibrium in the 
static de Sitter coordinates. Thus, the corresponding quantum state is determined by 
the Green function which is periodical along the time coordinate of the static frame with 
period if3. This enables one to use the Euclidean formulation of the theory. In particular, 
one can define the partition function Zp for arbitrary (3 as a functional integral, where 
the fields configurations belong to the spherical domain Si with conical singularities on 
the Euclidean horizon. 

There are at least two motivations for this work. The first one is to obtain the one- 
loop effective potential on such non-trivial spaces like St in an analytical form. This is a 
continuation of the pioneering studies by Shore |I[ and Allen [0, |3|] which investigated the 
effective actions on hyperspheres S 4 . Furthermore, the paper extends the results of Ref.[|J 
concerning the scalar fields, to more interesting gauge models with spontaneous symmetry 
breaking. A second motivation is related to the possible application of these results to 
cosmological problems since the exponentially expanding phase of the early universe |J is 
described by the de Sitter geometry. Our analysis has become possible after the spectrum 
of vector Laplacians on singular <i-spheres has been found explicitly in Ref. || . 

In this paper a special attention is paid to calculation and comparison of the effective 
potentials for two physically relevant quantum states: the de Sitter invariant (dSI) and 
the zero-temperature (ZT) states. They both can be considered as possible vacua in 
quantum field theory in the de Sitter space. 

The dSI state resembles the Poincare invariant vacuum of the Minkowski space-time 
since it preserves the symmetry of the de Sitter space. The corresponding temperature 
is equal to = (27ra)" 1 , where a denotes the de Sitter radius [0] and the Euclidean 

section of the de Sitter space is the hypersphere S 4 . In analogy with black hole physics 
this temperature can be called the Hawking temperature. In the Gibbons-Hawking path 
integral || the Euclidean gravitational action with a positive cosmological constant has 
the extremum on a 4-sphere S A . Thus, the choice of dSI state has a natural explanation 
in the semi-classical treatment of the quantum cosmology. For this reason the one-loop 
effective action for quantum fields on S A has been studied in a number of works |§-[[|. 
||, [10[| • m particular, some authors [[3], ||, [L0| investigated the phase structure of the GUT 



models with respect to the value of the de Sitter radius a. 

The dSI vacuum is analogous to the Hartle-Hawking vacuum ||TTJ introduced for quan- 
tum fields around an eternal black hole. From this point of view the zero-temperature 



state ((3 — > oo) in the de Sitter space is analogous to the Boulware vacuum |12| . The 



properties of this state are much less investigated, and we use here the explicit expressions 
of the effective potential to see how the choice of the vacuum affects the phase structure 
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of the gauge theories with spontaneous symmetry breaking. 

The paper is organized as follows. We briefly discuss in Section 2 why the Euclidean 
functional integral on the spherical domains Si can be interpreted as a partition function 
in the static de Sitter space. The FT effective potential for gauge models is computed 
with the ^-function method in Section 3. In Section 4 we compare the symmetry breaking 
mechanisms at different values of a in the two vacua for models of scalar electrodynamics. 
Summary of the results is given in Section 5. Technical details concerning the derivation 
of the (^-function are left for Appendix. 

2 Partition function 

In the static coordinates the line element of the de Sitter space with radius a can be 
written in the form 

ds 2 = cos 2 xdt 2 - a 2 (dx 2 + sin 2 xdOj + sin 2 x sin 2 Q\ddl) , (2.1) 

where — oo < t < +oo, \x\ < tt/2, and < 61,62 < tt. The above coordinates ( |2.1| ) 
cover only a part of the space. One can consider a Killing vector field generating a 
one-parameter group of isometries, subgroup of SO (1,4). Then coordinates ( |2.1| ) can be 
related to the time-like part of the Killing field associated with the translations along the 
time t. This region is restricted by the Killing horizon where the Killing field is null. It 
also coincides with the event horizon for the observers with coordinates x = 0. 

Let G(x, x') be a scalar Green function in the de Sitter space-time for dSI state (its 
some explicit expressions can be found in |13[-|[T5[). This function is periodic when the 
time coordinates of points x or x' are independently increased by i/3n where j3u = Ina. 
Thus the Euclidean Green function 

G E (T-r',u,u')=iG(i(T-T , ),u,u') , (2.2) 

which is obtained from G with the help of a Wick rotation, is defined on a 4-sphere S 4 
with the line element 

ds 2 = cos 2 xdr 2 + a 2 (dx 2 + sin 2 xd6 2 + sin 2 x sin 2 6 x d6\) , (2.3) 

where < r < j3n- Note that because of the periodicity property in imaginary time 
the function G(x, x') can be interpreted as a Green function for a canonical ensemble 
at temperature . The physical meaning behind this interpretation is that any freely 
moving observer experiences the dSI state as a thermal bath at temperature (3jj , see 
Ref.0. 

A natural generalization of G(x, x') is a Green function Gp(x, x') in the de Sitter space 
which is periodic in imaginary time with an arbitrary period j3. Functions of this kind 
can be constructed from G(x,x') with the help of a reperiodization formula suggested by 
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Dowker [I7|]. One can interpret Gg( Green function of a canonical ensemble 

of particles in the static part ( |2.1| ) of de Sitter space at temperature (3~ l . It is also possible 
to introduce the Euclidean function Gp(x, x') by using a definition analogous to Eq.( [2.2|) . 
The space where Gg(x,x') is set on is a spherical domain Sp, which is described by the 
metric ( |2.3|) of the 4-sphere, but now with period (3 along r. 

The main property of Sp is that it has conical singularities at the points \ = ±7r/2, 
near which the space looks as S 2 x C@, where Cp is a cone with deficit angle 27r(l — (3 / (5h)- 

Let us consider now the functional integral 

Z p = f [Z>0]e~ /j5[<M ^ ] , (2.4) 



where Ie[4>, 9^u] is the classical Euclidean action for the fields on St, g^ u is the metric 
tensor defined by Eq. (|2.3| ) and [Dtp] is a covariant integration measure. The function Gp 
can be obtained from this integral in the standard way. For instance, for a real scalar 
field <fi it reads 

Gf(x, x) = Zp 1 J [D<f)}<f)(x)<f){x')e- lE[(t> ' 9 ^ ] . (2.5) 

It follows from ( |2.4| ) and ( |2.5| ) that integral Zp is analogous to a statistical-mechanical 
partition function Tr e~@ H of a canonical ensemble of particles at temperature Thus 
we will call Zp partition function and use it for the definition of the effective potential in 
the theory. The corresponding quantum states will be called FT states. Some justification 
of this way of doing is that at /3 = /3# the integral Zp defines the quantum theory and 
effective potential in the dSI state. It should be noted, however, that in case of the Killing 
horizons Zp and statistical-mechanical partition function Tr e~^ H are not completely 
equivalent Q 

Their relation was investigated in the case of black hole geometries with the aim to 
establish the statistical explanation of the black hole entropy, see for instance Ref.s||18|- 



2*1| . We will not dwell on this issue further since it is not relevant for our purposes. 
It has to be also mentioned that in the presence of the Killing horizon the function 



Gp(x,x') does not have the Hadamard form [^] when (3 ^ (3r- In this case the stress 
energy tensor of a quantum field has a non-integrable divergence on the horizon of the 



chosen static coordinate system ||16|| . As was discussed in Ref . ||23|| , this property alone 
may be not sufficient to exclude FT states as unphysical, since the computation of the 
stress tensor neglects the backreaction effects. Such effects are very strong near the 
horizon and a reliable computation requires a nonperturbative approach which is not 
still developed. However, a non-Hadamard form of Gp(x,x') is not an obstacle for the 
definition of the partition function ( |2.4j) . The conical singularities of Sp result only to a 
number of additional ultraviolet divergent terms in the effective action — In Zp and one 
can use the standard renormalization procedure to give a meaning to the integral 
We will discuss this later. 



1 Note that this disagreement is also present at (3 — fin- 
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3 Effective potential 



3.1 The model 



We focus now on the definition and computation of the FT effective potential for the 
scalar electrodynamics in the de Sitter space. A similar computation has been done 
before by Shore JTJ and Allen |2j for the dSI state. We restrict our analysis to the abelian 
theories because the generalization to non-abelian case, as was shown in Ref. is almost 
straightforward. 

Let us consider the model of a complex self-interacting scalar field (f), which interacts 
also with an abelian gauge field A^. The classical Lorentzian action reads 

1 



I [A 



'— g dx 



where 



(Dyy - v{<f>) 



D, 



v{<t>) 



X 

+ 4! 



(3.1) 

(3.2) 
(3.3) 

(3.4) 



The coupling e stands for the electric charge, A is the constant of the self-interaction, which 
is assumed to be positive, and £ denotes the coupling to the scalar curvature R = 12 /a 2 . 



According to Eq. (|2.4j ), the partition function Zp for this model at an arbitrary temperature 
f3^ 1 is the Euclidean functional integral 



[D(j>] [DAJ [Dc] [Dc]exp{-/ £ L4 M ,0,c,c]} 



(3.5) 



I E [A^, ( f),c,c} = I E [A fM ,4>] + AI[A IM , ( f)}+I gh [c,c,^} , (3.6) 

where all the fields are given on the spherical domain Si, Eq.( [2.3| ). The functional Ie 
denotes the Euclidean form of action ( |3.1| ). The term A J is a gauge fixing term, c, c are 
the corresponding ghost fields and igjjc, c, 0] is their action. The explicit form of these 
quantities will be fixed below. 

We will use the effective potential method to study the phase transitions in the theory 
( |3.1| ). The symmetry breaking in our model is characterized by the average value < <p >p 
of the scalar field in the given FT state. If this average does not depend on time and 
on the spatial coordinates x, #2, the effective potential Vgg(v?, /3) can be introduced by 
following the method of Ref. |4j] . To this aim one can separate the field onto a constant 
part (p (which can be chosen real) and the excitations <fi' — <\>\ + i(f>2 



4> = tp + <p' = tp + 4>i + i<p2 ■ 

The effective potential is determined by the relation 

Z p = [dip exp{-/W eff (^,/?)} , 



(3.7) 



(3.8) 
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exp {-(3VV eS {<p, (3)} = J [D<f>>] [DA,] [Dc] [Dc] exp {-I E [A^ 0', c, c}) , (3.9) 

where (3V = ^na 3 j3 is the volume of Si. Note that the integral ( |3.8|) is the usual integral. 

As was shown in Ref.0, a point y?o where V e fl(tp, (3) has a minimum, coincides with 
the one-loop value of the average < (ft >p. The real part of V Q fi((p, (3) is a sum of the 
classical potential energy V(ip) and a quantum correction. If the field configuration <p is 
unstable, then V e Q((fo, (3) has a nonvanishing imaginary part which determines the decay 
probability of this configuration. 

To compute Vgg with the help of Eq . ( ^9|) it is suitable to use the t'Hooft gauge with 
the gauge fixing term 



-a 

2 J si 



Jgdx A \V ^ - aC l e V (j) 2 



(3.10) 



where a is an arbitrary parameter. Then, by expanding the classical potential to the 
second order in (pi and <fi 2 , one gets 

I E [A M <f>]+AI [A„ <I>] = \L[K (Ag? + eVV) A» T + c< (-V M V„ + ert W) ^ 

+<M~V 2 + V»)0i + 2 (-V 2 + a^eV + ^V'(0))0 2 ] ^d 4 x , (3.11) 

where Aj, and = V M x are t ne transverse and longitudinal component of the vector 
field, respectively, and Affi = —g fl ^'V CT 'V (T +R, L , is the transverse Hodge-deRham operator. 

Hereafter, we will be interested in the situation where the one-loop quantum effects can 
change significantly the form of the classical potential. As was pointed out by Coleman 
and Weinberg p^| , who investigated the model ( |3.1| ) in the flat space, one has to assume 
to this aim that the gauge coupling e 4 is of the order of A. In this case the quantum 
corrections due to the gauge fields will be comparable with the classical potential. On the 
other hand, to provide the convergence of the perturbation expansions one has to assume 
that A 1. In this case one can safely neglect the contributions to the effective potential 
of the scalar loops which will be proportional to A 2 e 4 . 



3.2 Representation for the potential 

In order to simplify the computations it is convenient to choose the Landau gauge for 
which a — > oo. Then the contribution of the ghost fields does not depend on the average 
value of the field tp, and since it does not affect the symmetry breaking in the model, it 
can be neglected. As was explained above, the contribution of the scalar excitations can 
be neglected as well. After that one arrives to the obvious result 

3 



V eS (<p,(3) = V(<p) 



log det 



^ 2 (Ag? + 



Further we make use of the ^-function regularization method [2 
terminant of an operator A with eigen-values A n 



(3.12) 
which defines the de- 



log det 



'-A 



-d(o)-Uo)io g // 



(3.13) 
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in terms of the generalized (^-function 



Ca(z) 



E K z 

n,(A»^0) 



(3.14) 



and its derivative C'a( z ) = dC,A{ z ) I 'dz. An arbitrary mass parameter \x is introduced in 
Eq.( |3.13|) to keep the right dimensionality. 

Let us note that the (^-function method automatically gives a finite expression for the 
quantum determinants. The last term in the right hand side of Eq.( |3.13| ) corresponds 
to the finite counterterms which are always present after renormalization. The conical 
singularities are known to introduce additional ultraviolet divergences, see for instance 
Ref's |2l| The renormalization on manifolds with conical singularities has been 

discussed recently in a number of papers |2R [2"8fl-f55fl,[jB]. 

In the case of the Hodge-deRham operator all one-loop divergences which are linear in 
the conical deficit angle, as shown in Ref.0, are removed under standard renormalization 
of the gravitational couplings []3JJ in the bare gravitational action. The renormalization 
of the remaining divergences of order (/3 — fin) 2 in the deficit angle or higher is more 
involved, and requires additional counterterms in the effective action. These terms have 
the form of integral invariants defined on the singular surface ||3U|| . The values of the 
additional couplings cannot be predicted in the theory, and so to simplify the analysis we 
will assume that these couplings in V q q are absent. 

The definition and representation of the ^-function ^ T \z) for the transverse Hodge- 
deRham operator is given in Appendix A. It is convenient to parametrize this function of 
z, in terms of the parameters (3 and a = j — e 2 a 2 ip 2 . Thus, Eq.( |3.12| ) reads 



]ira 3 (3 



U {T) (0,/?,a) + C (T) (0, (3, a) log 
dz 



(3.15) 



Our key result, whose derivation is given in the Appendix A, is that on the ^-function 
Q 7 ^ (z,0,cr) of the vector operator a 2 A^ + g^ u {\ — er) and the (^-function (z,j3,a) 
of the scalar operator — a 2 V 2 + 9/4 — 0" are related in a simple way, see Eq. ( |A.8|) . In 
particular, one can show with the help of Eqs.( A~9 ), ( A. 10 ) and ( A. 12 ) that for arbitrary 



d 

dz 



C (T) (0, 13, a) 



3|c<°>(o, M 



C (T) (0,/3,a) 



+ 
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u ) ip(u)du 



3C (0) (o, M + f§ 



a 



(3.16) 



(3.17) 



where i/j(u) = T'(u)/T(u) is the Digamma-function [55]. Thus, the quantum correction 
V^g to the potential due to gauge fields and the correction V£q from the scalars are related 
in the universal way 

V° s (a,f3) = 3V° s (a } f3)+{3- 1 n(a) , (3.18) 
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where the function Q(a) is temperature independent 



Ana 3 



+ 



— — it J ip(u)du 



1 

- I cr 

2 



11 \ , > > 

) kg/* a 



(3.19) 



Hence, the study of the effective potential ( pj.!5|) is reduced to the investigation of the 
scalar functional V^(a, /3) which has been done in Ref. Let us note, however, that 
even the structure of ^f( cr > 0) is rather complicated and in general one may rely only on 
numerical calculations. 

What is interesting is that V^(a, (3) can be found in an analytical form in the two 
most interesting limits: in the dSI state and in the ZT state. We make use of this fact 
to consider the phase transitions in the gauge model at different values of the de Sitter 
radius and compare the phase structures of the theory in these two cases. 



3.3 DeSitter-invariant state 



We assume that the renormalized mass of the field ip is zero and so the classical potential 
in Eq. ffTIgD is 

(3.20) 



with R = 12 /a 2 . Then the expression for the potential, which follows from Eq. ( [A. 13] ), is 



2 , ^..4 



(47r) 2 a 4 



1222,1444 19 

-e a ip + -e a ip 

3 y 4 y 192 



( -3, |) - §Ck 



"'2 



r— \ftr\ ( 3 s 

u I u — - ) (u — 3)4>(u)du 



+ (^ 4 aV + e 2 aV + ^) log A* 



(3.21) 



One can also define the energy E and the entropy S of the quantum fields at the Hawking 
temperature. Thus, by making use of ( |A.13|) , we find 



E 



Ana 3 d r 



\P=Ph 



Ana 3 ( 1 . _ 2 A 4 N 



Ana 



1 2 

-■-444 ^222 

-e a (p e a ip 

8 9 



(3.22) 



The entropy S of the quantum field in the dSI state is defined as S = 2na{E — V q q) and 
its expression follows from Eqs. ( [3.2 1|) and ( |3.22|) . 
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3.4 Zero-temperature state 

As it follows from Eg. ( |A.l7 ), the effective potential in this case reads 

v e s(<p, P = 00 



1 2 ^4 3 
2^ + 4! * - p¥ 



j- 4 4 4 



11 2 2 2 317 
— e a 09 H 

12 Y 192 



-3,|) + |c«(-l,|) + 3 



-\f5 



+ 



2 

-y/a 



1 



I it ( it — - ) [u — l)ip(u)du 



1I4441I222 1 
+ -e a r + -e a or 

\4 ^ 4 ^ 40 



u (u — 1) ip{u)du 

1 2 2 
log// a 



(3.23) 



It is interesting to note that the quantum correction to the potential due to gauge fields 
at zero temperature has the same form of the correction due to the scalar fields multiplied 
by factor 3, see Eq.( [3.18|) . 



3.5 Flat limit 

By taking into account that for z ^> 1 

'3 



t, , /3 \ , 11 

Re w — h iz\ = log z H z 

y V2 J 5 24 



Re ip 



iz 



log z 



127 
96C? 
7 



(3.24) 
(3.25) 



,2 / 24 960 
one can check that the energy density (47ia 3 /3)~ 1 E and the effective potentials V^g(v?, Ph) 
and V Q ff((p, (3 = 00) coincide in the limit a —>■ 00 with the effective potential in the 
Minkowski vacuum Vj^ink 



A 



3e 4 



%ink(^) = 4^ + 64^V' 



log 



(3.26) 



The last expression enables one to fix the unknown constant /x in terms of the physical 
parameters. 

The effective potential ( |3.26|) has a minimum at a non-zero value ip = <po where 

3 2 ' ^ 



2 Z 1 / 

^0 = ^2 exp 1 1 



9e 4 



(3.27) 



Thus the quantum state is characterized by a field configuration with the non-zero average 
< (f> >~ {p , where the gauge symmetry is spontaneously broken and the gauge field 
becomes a vector boson with the mass M = e\ipo\. The potential ( 3.26 ), written in terms 
of the parameter M, reads 



3e 4 



%inkM = g^aV 



log 



eV 
M 2 



(3.28) 



It coincides with the Coleman- Weinberg expression [p4j| . The result that in the flat limit 
dSI and ZT states coincide with the Minkowski vacuum is not surprising. In this limit the 
de Sitter group converts into the Poincare one and the Hawking temperature vanishes. 
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3.6 Limit of large curvatures 



It is also interesting to compare the form of potentials ( |3.21p and ( |3.23| ) in the opposite 
limit, when the curvature of the space-time is very large, namely e 2 a 2 (f 2 <C 1. Since we 
are interested in the phase transitions, terms in the potential which do not depend on 
if are irrelevant. Hence, from now on it is more convenient to deal with the difference 
V e ff((p,P) — V e Q(0, /3) which vanishes for ip = 0. 

By taking into account Eq.( |3.2lD one can easily prove that the dominant contribution 
for the large curvature at (3 = Ph is 



V eS (<p,p H )-V eS (0,(3 H ) 



3eV T n 19 

— P H 

167r 2 a 2 L 6 

3eV 4 



27 -log a 2 M 2 



647T 5 



3 -2 7 - log a 2 M 5 



where 



32tt 2 £ 8 



7T 



2 A 



9e 4 



3 
2 



(3.29) 



(3.30) 



and 7 = .5772.. is the Euler constant. As far as the ZT state is concerned, one can prove 
the analogous formula 



3eV 
Q4ir 2 a 2 
3eV 
64vr 2 



V eS (<p,P = 00) -V eS (0,P = oo) 



+ 



Q + fT-loga 2 M 2 

K + — - log a 2 M 2 
3 



where 



K = 12 



1/2 /l 







— z 







128vr 2 £ 8vr 2 A 14 
~~ ^ 9^ + ¥ 



cfe = -5.97011. 



(3.31) 



(3.32) 
(3.33) 



Equations (|3.29|) and (|3.31| ) show that in the limit when the de Sitter curvature is large 
the potentials are different but have similar structures. 



4 Symmetry breaking 



In order to compare phase transitions in the scalar electrodynamics (|3.1| ) in dSI and ZT 
states we rewrite expressions (|3.21|) and fl3.23j ) in the following dimensionless formQ 



Ai(x,y) 



64vr 2 
3M 1 



v eS &,p H )-v eS (o,(3 H ) 



4- P + — -logy 
y V 6 



2 ! J_ 4 

x log y + — 

y 



u [u — — ) (u — 2>)ip{u)du 



(4.1) 



2 For dSI state representation (4.1) was suggested in || 
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Mix,y) 



64tH 
3M 1 



V eS fo>, (3 = oo) - F eff (0, = oc) 



x 



-(Q - logy) + x 2 

y 




+ 



uxu-- 



[u — \)ip{u)du 



u(u — l)ip{u)du 



Here 



x 



eV 
M 2 



y = a 2 M 2 



a 



4~ Xy 



(4.2) 



(4.3) 



and the renormalization parameter /x has been expressed, according to Eq. (|3.27|) , in terms 
of the vector boson mass as fi 2 = M 2 exp — lj- Interestingly, all the information 
about the constants e 2 , A and £ of the considered model is contained in the parameters 
P and Q defined in Eqs.( p.30| ) and (|3.32|) , respectively. 

For arbitrary values of constants e 2 , A and £ the parameters P and Q are independent. 
They can be related to each other only in the case of the minimal coupling £ = 0. 

We consider now the phase structure of the theory in the limit of large and small 
curvatures of the de Sitter space, where it can be investigated analytically. 



4.1 Phase transitions at small curvatures 



In this limit e a ip = xy 3> 1 and one has 



At(x,y) = x 2 ( 



log x 

& 2 



4- logx + P+ - 
y V 2 



A 2 (x,y) = x 2 



logs - ^ + - (\ogx + Q - 



(4.4) 
(4.5) 



X 

lloga; + g- 

y V 3 

Eq. (|4.4f) coincides with the result of Ref. ||. As was shown in Ref. the symmetry 
breaking in this limit in dSI state is characterized by first order phase transitions. Under 
such transitions, when the curvature radius a reduces to a critical value a c , the mass of 
the vector boson changes discontinuously from a value M c = eip c to zero. By comparing 
Eqs. ( J4.4Q and ( [4.5|) one can see that the effective potential at zero temperature has a 
similar structure, so in this case one may expect an analogous behavior. 

Let x c and y c be critical values of parameters x and y in the point of the phase 
transition. They are related to the critical radius a c and mass M c as follows 

M 2 

^=^| , Vc = a\M 2 (4.6) 

where, as before, M is the mass of the vector boson in the Minkowski space-time. In the 
given approximation the values of x c and y c , with x c > 0, can be found explicitly. They 
must satisfy the following conditions 

d 



dx 



A(x,y) 







(4.7) 



x=x c ,y=y c 
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which can be resolved. For dSI state the results coincide with those of Ref. 



Mj = M 2 exp 



2M: 



.5 _p 
2 1 



1 + P + \/P 2 - 4 



P + VP 2 - 4 



For the zero temperature we find the following critical mass M c and radius a c 



M 2 C = M 2 exp 



f-Q 



Q-f + \l(Q-f) 



~a 2 . 



-M~ 2 
2 c 



Q 




(4.8) 

(4.9) 
(4.10) 

(4.11) 
(4.12) 



Eqs. O and (|4TT0|) are valid for P > 2, and Eqs. ( pi]) and ( [lip for Q > 49/6. 
Thus we can conclude that in the limit of small curvatures and for the given values of 
P and Q the dSI and ZT states have qualitatively similar properties. In both cases one 
has first order phase transitions, but with different values of critical radii and masses. 
These properties are shown in Figures 1 and 2, where the functions Ai(x,y) and A 2 (x,y) 
are evaluated for P = 10 and Q = 16, respectively. The numerical computations of the 
critical parameters are in agreement with expressions (|4.9| )-( [lT2| ). 



4.2 Phase transitions at large curvatures 



As one can see from these 



We can use in this regime asymptotics ( |3.29|) and (|3.31[ ! 
equations, the terms proportional to ip 2 change the sign from positive to negative when 
a becomes larger than a critical value a c . In this case to have the symmetry breaking 
in the theory at a > a c the terms ip 4 in (|3.29|) and (|3.31|) have to be positive. When a 
becomes smaller than a c the mass of the vector boson gradually vanishes. This situation 
corresponds to second order phase transitions. 
The critical radius a c in dSI state is @ 
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Mr 



cxp 



r ^ is 



(4.13) 



In order to have a second order phase transition one must ensure that the coefficient of 
the <^9 4 term in ( |3.29|) is positive. By substituting (^4.13|) in Eq.( |3.29|) we find that the 
above condition is satisfied || for 



3-2 T - log a\M 



-P 



6 



> 



P < 



1 

6 



(4.14) 



The value P = P„ 



•1/6 represents in fact a critical boundary which separates the 



models of the scalar electrodynamics with different kinds of phase transition. For models 
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with P above P cr the system undergoes first order phase transitions, whereas for models 
with P below P cr one has second order phase transitions. The value of P cr is also confirmed 
by the numerical analysis. 

A similar expression for the critical radius a c at zero temperature can be found from 
Eq.flOID 

a\ = M~ 2 exp {Q + K} . (4.15) 
The difference between the values of the critical radii can be estimated by their ratio 

S =etp (^ + Jr + 3 + 2 7+ ^!) . (4.16) 
a 2 \ e 2 e 2 J 

It depends only on £ and e 2 , and for £ ^> e 2 one can conclude that a c 3> a c . For the 
minimal coupling (£ = 0) a? c /a 2 c ~ 1.63. 

In the case of ZT state the critical value Q cr can be determined as well. By substituting 
Eq. (|4.15|) in the coefficient of the quartic term of Eq. (|3.31| ) one gets 



2D 2D 2D 

K + j -loga 2 M 2 = |-g>0^g<| . (4.17) 

Thus Q cr = 20/3 and it is confirmed by the numerical computations. 

The typical behavior of the effective potentials in case of second order phase transitions 
is illustrated by Figures 3 and 4, which are obtained by using formulas (|4.1|) and (|4.2|) . 
Figure 3 depicts the function A\(x, y) for P = — 1/2 and Figure 4 gives A 2 (x, y) for Q = 6. 
In both cases the critical radii a c and a c are in agreement with Eqs. (|4.13|) and ( |4.15|) . 

4.3 Classification of models 

As we have shown, the first order phase transitions always take place in the models of 
scalar electrodynamics (|3.1| ) when constants P and Q are large and positive. On the other 
hand, if these constants are negative with large absolute values, the symmetry breaking 
corresponds to second order phase transitions. 

An analytical study of the region of parameters £, A and e 2 (or P cr and Q cr ), where 
the kind of the phase transition changes from the first to the second order is difficult 
because the effective potentials have quite involved forms (|4.1|), ([472]). Thus, in principle, 
one should use here numerical methods. It is interesting, however, that the analytical 
estimates P cr = — | and Q cr — y of Section 4.2 are in very good agreement with the 
numerical results. For dSI state this was first pointed out by Allen 0. 

As we already mentioned, the symmetry breaking mechanism in the considered models 
is completely determined by the parameters P and Q. By using their definitions ( |3.30|) . 
( |3.32D one finds that 

«^+^+! . 

Note that the parameter A must be chosen positive in order to have a scalar potential 
bounded from below. We also assumed that A ~ e 4 in order to have a considerable 
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quantum correction to the classical potential coming from the gauge fields. For values 
A ^> e 4 our perturbative approach is not reliable, so we have to restrict ourselves to the 
interval of parameters < Ae~ 2 < e 2 . Thus, by taking into account that e 2 /4n ~ 1CT 2 , we 
find with the help of Eq.(g7L8]) the allowed region in the P-Q plane < Q-4P-32/3 < 1. 
There we can make definite predictions about the properties of the considered models. 
This region is shown on Figure 5. It consists of three subregions (I), (II) and (III), bounded 
by dashed lines. Each subregion indicates a class of models with the same properties. 
Region (I) is determined by inequalities P < —1/6 and Q < 20/3. If the coupling 
constants £, e, A satisfy these restrictions, the corresponding models undergo second order 
phase transitions in both dSI and ZT states. 

In region (II) one has P > —1/6 and Q > 20/3 and first order phase transitions in both 
dSI and ZT vacua. 

Region (III) is the most interesting, since in this case P < —1/6 but Q > 20/3, thus the 
system undergoes a first order phase transition in ZT case, whereas it has a second order 
phase transition in dSI case. 

In case of minimal coupling (£ = 0), according to Eqs. (|3.30| ) and ( ft.32| ), P < — 3/2 < 
—1/6 and Q < 14/3 < 20/3, and so only second order phase transitions are possible in 
both states. 



5 Summary 

We investigated the effective potential for the gauge models on the singular spherical 
domains Si. The spectrum of the corresponding wave operators on these backgrounds 
can be found exactly and it enables one to calculate the potential with the help of the 
^-function. We have shown that for arbitrary values of f3 there is a simple relation ( |3.18| ) 
between the one-loop corrections to the potential from gauge and scalar fields. Thus in 
studying the gauge models in the de Sitter space one can make use of the computations 
of the scalar potential of Ref.jp. The effective action on Si can be interpreted as the free 
energy of fields in static de Sitter space, where the parameter (5 is the inverse temperature 
of the system. Thus, by changing (3 one can learn how the presence of temperature affects 
the properties of such a quantum theory. The aim of the paper was to compare the phase 
structures of the gauge models in the de Sitter invariant vacuum, which is usually used 
in the quantum cosmology, with the structure of the normal vacuum which is the state 
at zero temperature. 

According to our analysis both vacua have very similar properties in the extreme 
regimes of very small and very large curvatures. If the restoration of the gauge symmetry 
happens at large radius a (small curvature) of the de Sitter space, then the corresponding 
phase transition is of first order regardless what vacuum is considered. Analogously, if the 
symmetry is restored at small values of a (large curvature), one has second order phase 
transitions. The critical masses and radii are completely determined by the constants P 
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and Q, which are the functions of the parameters £, e 2 and A. Expressions (|4.11|) , ( |4.12| ) 
and ( |4. 15| ) for the critical parameters in ZT state represent the new result. 

We have also found the critical value Q cr = 20/3 for the parameter Q , which separates 
models with second (Q < 20/3) and first (Q > 20/3) order phase transitions in ZT state. 
Finally we proved the existence of a class of the models (P < —1/6 and Q > 20/3) where 
the kinds of phase transition in the two vacua are different. 

The detailed analysis of the effective potential for all /3's is not given here, but it can 
be carried out in principle by making use of the properties of the ^-function. In the most 
interesting cases, however, it is sufficient to restrict the study to the expansions near the 
zero and Hawking temperatures, as it was done in Ref.jp. 

Acknowledgements: This work was supported in part by the Natural Sciences and 
Engineering Research Council of Canada. 
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A (^—function 

Here we consider ^-function 



C (T) (z) = EE^(Ay^ + E D T n (A T nfi ) 

n=0 1=1 n=l 



(a.i; 



for the transverse operator in Eq.( p.l2[ ) which is defined for convenience in the dimen- 
sionless form as a 2 (A^ + g^e 2 ^ 2 ). 

According to M the degeneracies in four dimensions read 



n 



Ztf=3(n + l)(n + 2) , D 1 n = -(3n + 5) 



(A.2) 



while the eigen-values A^ are 

= (n + ll) {n + jl + 3) + 2 + e 2 a V , (A.3) 

where 7 = Ph/P- As in Section 3, it is suitable to parametrize the function ( |A.1|) by the 

1 
1 



parameters P and a — j — e 2 a 2 (p 2 . Thus one can write 



71=0 1 = 1 



n=l 



where 



AL(a) = (n + 1 l + 3/2) 2 -a . 



(A.4) 



(A.5) 



Function ( |A.4| ) has an interesting relation to the ^-function of the scalar Laplace operator 
— a 2 V M V,j + 9/4 — a on the same spherical domain Si. These operators were studied in 
fl, where they were shown to possess the eigenvalues coinciding with ( |A.5| ). So the scalar 
^-function is represented in the form 



(^(z,P,a) = EE^ 0) (Al,(a))- 2 + E^ 0) (Al,o(-))- 



(A.6) 



n=Q 1=1 



n=l 



where the degeneracies 



£>W = (n + l)(n + 2) , 5i 0) = \{n + l)(n + 2) 



(A.7) 



differ from the degeneracies (|A.2|) of the vector operator. By comparing Eqs. (|A.6| ) and 
( |A~7D with (|AT|) and ( |A~2|) one obtains 



C (T) (z,P,a) = 3((°\z,P,a)-2Y: 



n=0 



n H — 
2 



3 

n H — 
2 



<7 



3C {0) (^, M-2£ ^C,(z)Cfl 2z + 2* - 1, - 

fc=0 v 2 



(Ai 
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where Cy.(z) = T(z + k)/(k\T(z)). To get the second equality of this equation one has to 
decompose the last term in a and make use of the definition of the Riemann zeta-function 
Cr{z, a). The function C (2, °") h as been studied in detail in the previous works [Q. [|. 
and this makes the analysis of Q^ T \z, (3, a) much more easy. In particular, by taking into 
account the result of [§J and Eq.( |A.8|) one obtains 



C (T) (0,(3,a) = 3(^ (0,(3,a) + ^--a 



11 



12 



3 

7 



51 - 6O7 2 - 87 4 2 7 



2880 



3 a 2 

—a H 

24 12 



11 
12 



a 



However, the analysis of the derivative of the generalized (^-function 

U {T) (0,(3, a) = 3^ (0,(3, a) 
dz dz 



- 2 



2^-1 



(Tip 



00 k / o 



(A.9) 



(A.10) 



is more involved and we consider only two of the most interesting cases where the result 
can be obtained in the simple analytical form. 

- Temperatures close to the Hawking temperature - 

This case was investigated by B. Allen in Ref. and our aim here is to rederive his 
result for j- z C, {0) (0, (3 H ,a) in a different way with the help of Eq. ([A.8|) . For the scalar 
^-function one can find the following representation [[| 



(u — l)ip(u)du 



a a 
Y2 + T2 



+ 



Cr -3, 



4^ 



'2 



where ip(u) = T'(u)/T(u). By using relation Cf?( n + l)- 2 ) 
possible to show that 



(A.ll) 

iy^£,^(z)/n\ it is 



/3\ » a k+l 



( R [2k + 1,- 



--u)ip(u)du . (A.12) 



Thus by substituting Eqs. ( |A. and ( |A. 12|) into (|A.8|) and using the properties of the 
■0-function one gets 



dz 



C (T) (0,(3, a) 



-\f5 



( 3 \ / o \ / / \ 1 O- 2 11 

u\u — — \\u — 3)ip{u)du H — : — a 



+2Ck ( -3, ~) - \c'n 



41 _ o*_ _ 
144^ ~ ~8 ~ 5760 



4 2 i 

973 
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h ( 16a 2 - 40a + 9) U (| + y^) + ^ (| - ^))] + O ((/3 - (3 H f) . (A.13) 



192 v / V V2 v / V2 

As one can see, this expression coincides at (3 — @h with the result reported in Ref. 



- Vanishing temperature - 

To derive d^^/dz in the zero temperature limit one can use the method suggested in |4j]. 
Unfortunately, the final formula given in Q has a wrong form because of a misprint, so 
here we use the opportunity to represent the correct answer. According to relation (K 
and Eq.(B2) of Ref.g, one has 

P H lim ((3- 1 ^-C (T) (0,(3,a))=3(3 H \\m (/T 1 #C (0) (0, 0, a)) 



3±C k (z) { ^ )2k 

k=0 



dz 



dz 



2z + 2k-l 



( R (2z + 2k-3, 3/2) - -( R (2z + 2k-l, 3/2) 



3f(z,a) . 
(A.14) 

The function f(z, a) was introduced in [|], where it was shown to be related to {z, (3h, c), 
see ( [A.ll] ), by the differential equation 



ds/a 



f(z,a)(^) 2z - 1 ]=3cr z -\^(0,(3 H ,a) . 



(A.15) 



Thus, one can write 



d_ 

dz' 



/(0, a) = yfd ^ ^ [ 3 (^C (0) (0, (3 Hi V 2 ) - ^C {0) (0, (3 H , 0)) - 2(/(0, y 2 ) - /(0, 0)) 



3* C (o) (0,^,0) -2/(0,0) 
. dz 



(A.16) 



By making use of ( |A.11[ ) in ( A.16|) one finally finds 



ft, Urn (/3 _1 ^C (T) (0,/9,<r) 



3 1/' +/' " U^-^l(ii 



l)ij){u)du — 3V0- 1 / — / J — l)ij){u)du 



12 8 160 Si? V 2) 2 S H '2 



(A.17) 
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Fig. 1 
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Figure 1: The above curves, from the higher to the lower, represent the function Ai(x,y) 
for P = 10 and jy = 7.4, 7.8, 8.49, 9 and 9.5, respectively (see Ref. @). The critical 
value of de Sitter radius (y/y ~ 8.49) corresponds to the first order phase transition. 



Fig. 2 
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Figure 2: The curves, from the higher to the lower, illustrate A 2 (x, y) for Q = 16 and 
■Jy = 3.5, 3.8, 4.2, 4.7 and 5, respectively. The first order phase transition takes place at 
^-4.2. 
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Fig. 3 



0.2 



0.15 



0.1 



0.05 



-0.05 



-0.1 




Figure 3: The curves, from left to right, represent the plots of Ai(x,y) for P = — 1/2 and 
y/y = 1.95, 2.13, 2.3, 2.5 and 2.8, respectively, (see Ref.0). The critical value yjy ~ 2.13 
corresponds to the second order transition. 



Fig. 4 
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Figure 4: The curves, from left to right, shows the function A 2 (x,y) for Q = 6 and 
Jy = .9, 1, 1.1, 1.3, 1.5, respectively. The second order phase transition happens at 
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P 



Figure 5: The regions (I), (II) and (III), bounded by dashed lines, represent the allowed 
areas in the P—Q plane. If the parameters P and Q are chosen in regions (I) or (II), one 
has always second or first order phase transition, respectively, in both dSI and ZT states. 
In region (III) first order phase transitions take place in ZT state, whilst dSI state shows 
second order phase transitions. 
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